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Tidal Disruption of Satellite Galaxies in a Semi- Analytic 
Model of Galaxy Formation 



o 
o 

Oh 

c/3 



O 

U 

Oh 
6 

o3 



> 

O 

in 

(N 

ON 

o 

On 
O 



Bruno M. B. Henriques 1 ' 2 *, Peter A. Thomas 2 

1 Institute of Cosmology and Gravitation, University of Portsmouth, Portsmouth POl 3FX, United Kingdom 

2 Astronomy Centre, University of Sussex, Falmer, Brighton BN1 9QH, United Kingdom 



Submitted to MNRAS 



ABSTRACT 

We introduce a new physical recipe into the De Lucia and Blaizot version of the 
Munich semi-analytic model built upon the Millennium dark matter simulation: the 
tidal stripping of stellar material from satellite galaxies during mergers. 

To test the significance of the new physical process we apply a Monte Carlo 
Markov Chain parameter estimation technique constraining the model with the K- 
band luminosity function, B — V colours and the black hole-bulge mass relation. The 
differences in parameter correlations, and in the allowed regions in likelihood space, 
reveal the impact of the new physics on the basic ingredients of the model, such as 
the star-formation laws, feedback recipes and the black hole growth model. 

With satellite disruption in place, we get a model likelihood four times higher 
than in the original model, indicating that the new process seems to be favoured by 
observations. This is achieved mainly due to a reduction in black hole growth that 
produces a better agreement between the properties of central black holes and host 
galaxies. Compared to the best-fit model without disruption, the new model removes 
the excess of dwarf galaxies in the original recipe with a more modest supernova 
heating. 

The new model is now consistent with the three observational data sets used to 
constrain it, while significantly improving the agreement with observations for the 
distribution of metals in stars. Moreover, the model now has predictions for the intra- 
cluster light, a very significant component of large groups and clusters, that agree with 
observational estimates. 
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1 INTRODUCTION 

In recent years, semi-analytic (SA) models have experi- 
enced a significant degree of success, achieving a consider- 
able agreement with a large set of observational properties. 
These range from the luminosity to the stellar mass func- 
tions, Tully-Fisher relations, clusterin g, galaxy colours, etc 
l|Bower et 'alj|200d ICroton et alj|2006h . 

Some fundamental changes determined the improve- 
ments achieved. First of all, these recipes were imple- 
mented on top of a direct dark matter simulation, which 
just recently have been able to simulate cosmological vol- 
umes with a large enough resolution to foll ow haloes con- 
tainin g dwarf galaxies (the Millennium run, ISpringel et al.l 
2005). Secondly, some key physical recipes were introduced: 
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the energy feedback from supernovae explosions, which 
for small enough galaxies and strong enough explosions 
can drive the gas o ut of the galaxy (jBenson et all 120031 : 
|Pe Lucia et al.ll2004l ); and the feedback from central black 
holes 



ies (IGranato et al.l 


2004 iBower et al. 1 


2006|; 


Cattaneo et al. 


2006 


; Croton et al. 


20061: Menci et al. 


2006 


: 1 Monaco et al. 


2007 


: Somerville et alj|2008h. 



The level of agreement achieved opens up a number of 
possibilities. This agreement can now be quantified against 
different observational data sets using robust stati stical tests 
l|Kampakoglou et al. I l2008l : iHenriques et al. 2009, hereafter 
Paper 1). This means that the allowed likelihood regions in 
parameter space can be obtained and confidence limits for 
the preferred parameter values can be built. Moreover, with 
the physics that determine the global properties of galaxies 
reasonably well understood, modellers can now focus on ad- 
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ditional ingredients that determine the fine tuning of galax- 
ies properties. 

Amongst these detailed st udies we can find the stripping 
of gas from satellite galaxies (IF ont et al.ll2008h ; the impact 
of the assumed dust mode l on the overall galaxy properties 
l|De Lucia fc Blaizot i20ojh and on the galaxy clustering over 
redshift (jGuo fc White! 120091 ) ; the impact of quasar mode 
feedback on the luminosity -temperature (L-T) relation in 
clusters ijBower et al.ll200a ): and the effects of a dynamical 
treat ment of galactic winds on host galaxies (|Bertone et al.l 
12007ft . 

Another physical process that is becoming more rel- 
evant in theoretical models is the disruption of stel- 
lar material from merging satellites due to tidal strip- 



semi-analytic 


dBullock et all 2001 


; iTavlor & Babull l200ll; 


Benson et al. 


2002; iMonaco et al 


20071; iHenriaues et al.1 


20081; ISomerville et al. 2008; Seek Kim et alj|2009l) and the 



l2009l ; IWetzel fc Whitdl2009h . 

Evidence for the importance of this process in galaxy 
formation comes from various fields. The existence of a 
dif fuse population of intra-cluster stars was first proposed 
bv IZwickvl rtl95lh a nd has since been detected unambigu- 
ously (iDurrell et aljlifjoi ; iGal-Yam et~aH 120031 : iNeill et all 
120051 ; iKrick et all 120061 ). The light associated with intra- 
cluster stars, or diffuse intra-cluster light (ICL), can con- 
tribute between 10 and 40 per cent of the optical emission 
of rich galaxy groups and clusters ([Bernstein et al.l [l995l; 
Gonzalez et al.ll2000l; IFeldmeier et al.ll2002l; IFeldmeier et al.l 
20041 ; iGonzalezet al 1120051 ; IZibetti et al.ll2005l ). 

Rather than having formed in the intra-cluster medium, 
gas-dynamical simulations generally agree that the bulk 
of the ICL is emitted by stars that have been continu- 
ally stripped from member galaxies throughout the life- 
time of a cluster, or have been eject ed into intergalac- 
tic space by merging galaxy groups ([Moore et al.l 19961; 



Napolitano et al. 2003; 



iMurante et al 
— T — I I — 
120051; 



|2004|; IWillman et all 
120061 ; 



20041; Sommcr-Lars en et al.l 120051; IMonaco et al.l 
Murante et al.ll2007l ; lRudick et alj|200gh . 

Observationally, low surface brightness features 
have been identified in the Coma and Centaurus 



clust ers dGregg fc West! 1 19981 ; iTrentham fc Mobasherl 
1 19981 ; IFeldmeier et all I2002D . indicating the pres- 
ence of dynamically-young tidal structures produced 



by the disruption of 



infalli ng galaxies. Moreover, 
(2005) show that the pro- 



iFaltenbacher fe Mathews! 
jected number density profile of dwarf galaxies in NGC 
5044 can only be explained by assuming that a significant 
amount of mass in satellite galaxies is tidally disrupted. 

If, as the evidence suggests, intra-cluster stars are the 
remnants of material stripped from merging satellites, the 
build up of this component can be directly followed by a 
semi-analytic mo del with a self consistent implementation of 
tidal disruption dBullock et al.ll20qil; ITavlor fc Babulll200ll; 
Benson et al] |2002| ; IMonaco et al.l 120061 ; ISomerville et al.1 



2008). Also, the disruption of material from satellite galax- 



ies might play a crucial role in solving a common problem 
to most semi-analytic models: the excess of d warf galaxies 
l|Henriques et al 1 l2008l ; IWeinmann et al.ll2006l) . 

Despite the fact that a number of authors have intro- 
duced satellite disruption in theoretical models of galaxy 



formation, the complexity of these models makes it difficult 
to fully understand the impact of different physics on the 
global galaxy properties. Also, the number of free parame- 
ters makes it some times impossible to take full advantage 
of the new physics introduced, when a manual tuning fails 
to reach the best possible fit. 

In this paper we make use of the Monte Carlo 
Markov Chain (MCMC ) statistical techniques introduced in 
IHenriaues et al.l l|2009l ) to test the impact of satellite dis- 
ruption in semi-analytic models of galaxy formation. With 
the introduction of a new physical recipe, both the normal- 
isation and shapes of the acceptable likelihood regions con- 
tain useful information that can be used to discriminate be- 
tween models. We will show that the new model is favored 
by the data, that it provides a better fit to the metallicity 
of galaxies than our previous model, and that it predicts 
ICL fractions that agree with observations. We note that 
similar MCMC tech n iques were proposed independently by 
iKampakoglou et al.l (|2008l ) to test the impact of different 
star formation modes in their specific semi-analytic recipe. 

This paper is organized as follows. In Section 2, we de- 
scribe the original semi-analytic model used in our study 
and the new implementation of satellite disruption. Section 
3 presents results from a model with disruption using the 
same parameters as in the original recipe to directly test 
the differences between the two. In Section 4 we briefly de- 
scribe the MCMC techniques introduced in Paper 1. In Sec- 
tion 5 we present the best fit model with satellite disruption 
and compared its likelihood maximum, acceptable likelihood 
regions and galaxy predictions with the previous model. Fi- 
nally in Section 6 we summarize our conclusions. 



2 THE MODEL 

To perform the work described i n this paper we used the 
Munich SA recipe as desc ribed in lCroton et al.l (120061 ) and 
|Pe Lucia fc BlaizotJ l|2007l . hereafter DLB07). It is built on 
top of a direct numerical simulations of the dark matter 
struc ture in a cosmologica l volume, the Millennium Simula- 
tion jSpringel et al. 2005). This simulation traces the evo- 
lution of dark matter haloes in a cubic box of 500ft" ^pc 
on a side. It assumes a ACDM cosmology with parameters 
fi m = 0.25, fi b = 0.045, h = 0.73, Q A = 0.75, n = 1, and 
as = 0.9, where the Hubble parameter is Ho = 100 /i _1 km 
s _1 Mpc -1 . The simulation follows 2160 3 dark matter parti- 
cles of mass 8.6 x 10 s /i _1 Mq. Since dark matter haloes are 
required to contain at least 20 particles, the minimum halo 
mass is 1.7 x 10 10 /i _1 Mq, with a corresponding baryonic 
mass of about 3.1 x 10 9 /i~ 1 M . 

The galaxy formation model itself follows the evolution 
of baryons from when they collapse into a hot gas phase, 
through cooling onto a disk where stars can form. As the 
most massive stars die, supernovae (SN) eject energy into 
the surrounding medium, reheating the cold gas back into 
the hot phase or even ejecting it into the external reservoir. 
The black hole evolution is modeled and in massive galaxies 
the mechanical heating from its quiescent growth suppresses 
the cooling. Mergers generate star formation bursts and, de- 
pending on the mass ratio between the galaxies, disks are 
destroyed to form bulges. Finally dust and stellar popula- 
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tion synthesis models transform the predicted quantities into 
galaxy properties that can be compared with observations. 

With all the recipes in place, the model is able to pre- 
dict reasonably well the shape of the luminosity function 
in different bands, the stellar mass function, the galaxy 
colours, the black hole-bulge mass relation, the metallicity 
of stars, the Tully-Fisher relation, etc. However some chal- 
lenges remain. For example, as happens with most current 
SA recipes, this model has known proble ms in correctly pr e- 



IWeinmann et al1l2006l ; iHenriques et al.l 



i y 

dieting the properties o f dwarf galaxies ( Croton et al.l 120061 ; 



20081 . 120091 ) and it 



has no predictions for the intra-cluster light, a component 
that can be substantial for the largest groups in the local 
Universe. 

These two aspects lead us to introduce a new physical 
recipe: the stripping of stellar material from satellite galax- 
ies due to the tidal forces they experience from the central 
companions. By removing this material from merging galax- 
ies we expect to decrease their size and hence the number 
density of low-mass galaxies (where an excess is seen in the 
original model). At the same time, by moving the disrupted 
material into the intra-cluster medium (ICM) we will be 
able to predict the properties of the 1CL. Moreover, since 
this corresponds to stars that would otherwise end up in the 
central galaxies of large groups and clusters, we expect to see 
considerable changes in the properties of the most massive 
objects in the model. 

To intro duce the new phy si cs we use the implementation 
proposed by iTavlor fe Babul l|200ll ) and follow the radial 
position of satellites as they merge into central objects. At 
each location, we compute the radius relative to the satellite 
centre at which the sum of the tidal force from the parent 
halo and the centrifugal force equal the self-gravitational 
force of the satellite. Any material outside this radius at 
each time step is considered to be unbound and is moved 
from the satellite galaxy into the ICM. 

We assume that disruption of stellar material only acts 
on galaxies that have experienced large enough tidal forces 
to completely strip their dark matter component (for those 
who are aware of the Munich SA terminology, type 2 galax- 
ies). 

2.1 Dynamical Friction 

The original model of DLB07 follows the merging of satel- 
lites by simply setting a up a "merging clock", whenever a 
satellite loses its dark matter halo. After this time is elapsed 
the galaxy is assumed to merge with its central companion. 

The mer ging time is calculated using the 

Chandrasckhar (1943) formula for the dynamical 
friction force acting on the satellite, as described in 



(1) 



iBinnev fc Tremaind (fl987h 



47rG 2 mL m(A)pB(x) 



where m sat represents the satellite mass, In (A) is the 
Coulomb logarithm (assumed to be ln(l + mh a io/w sat )), p 
is the local density and v aa ,t represents the orbital velocity 
of the satellite within the halo. B(x) is the error function 
and we take |/\/2crhaio = 1, where <7 ha io represents 

the 1-D velocity dispersion of the halo. 

Assuming circular orbits and an isothermal distribution 



of the total mass in haloes, the previous equation can be 
integrated to give the merging time: 



tdf ~ 1.17 



Gm sat ln(A) ; 



(2) 



where r sat is the halocentric radius of the satellite at the 
time that it lost its dark matter halo and V v \ v is the virial 
velocity of the central halo. 

For our implementation we use the same equation, but 
we now record the orbital radius of the satellite as it falls 
towards the halo centre. In this way, at each time step, we 
can compute the forces acting on the satellite at that radius 
and determine the amount of stellar material that becomes 
unbound. 

We note that the drag force acting on satellites and 
causing them to merge with central objects depends on their 
mass. Because of this, despite us using the same equation for 
the dynamical friction as in DLB07, the implementation of 
satellite mass loss will change the predicted merging times. 
As satellites lose their stellar material the drag force they 
experience will decrease, causing them to survive longer be- 
fore being accreted by the central galaxy. 



2.2 Tidal Disruption 

The mass loss on satellites occurs through the action of tidal 
forces. Assuming a slowly- varying system (a satellite in a cir- 
cular orbit) with a spherically-symmetric mass distribution, 
material outside the tidal radius will be stripped from the 
satellite. This radius can be identified as the distance from 
the sate llite centre at which the radial forces ac ting on it can- 
cel out l|Kindl 19621 : IBinnev fc Tremainelll987h . These forces 
are the gravitational binding force of the satellite, the tidal 
for ce from the c entral halo and the centrif ugal force. Using 
the King (1962) model and as described in lTavlor fc Babul 
(2001) the disruption radius r t will be given by: 



n 



Gm s . 



u4 t - d^/dr* 



1/3 



(3) 



where o; sat is the orbital angular velocity of the satellite and 
<f> represents the potential of the halo. 

Using the isothermal sphere approximation for the mass 
distribution of the central halo and satellite galaxy, and as- 
suming that the satellite follows a circular orbit, equation [3] 
can be rewritten as: 



n 



1 Caat 



y/2 cr ha i 



f s a t , 



(4) 



where <T sat and o"h a io are the velocity dispersions of the satel- 
lite and the halo, respectively, the former being estimated 
just before it becomes stripped. 

The material outside this radius is assumed to be dis- 
rupted and becomes part of the ICM. We assume that the 
galaxy as a uniform metallicity distribution so that equal 
fractions of stellar mass and metals are stripped away from 
the galaxy. We keep following the satellites until their or- 
bital radius becomes smaller than the sum of the radii of 
the central and satellite galaxies, at which point we assume 
a merger has occurred. 

To compute the forces acting on satellite galaxies, we 
assume that the total mass follows a spherically-symmetric 
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isothermal sphere. However, to determine the amount of ma- 
terial outside the disruption radius that is lost at each time 
step, we use appropriate distribution models for the stellar 
mass in the different galaxy components. The model nat- 
urally defines galaxies as a combination of a disc and a 
spheroid. We now explain how these two components are 
modelled. 



2.2.1 Mass Distribution in Disks 

We model the disc component as an exponential mass dis- 
tribution given by: 

rR 

M disk (<i?)=/ 2nRY. dfi e X p- R/Rd dR, (5) 
Jo 

where Ed,o gives the surface density of the disk and the disk 
sc ale- length ( Rd) is calculated using the formalism derived 
by I Mo et all (1 19981 ) . Integrating equation [S] the disk mass 
inside a given radius R is given by: 



M disk (< R) = M disk 



R_ 

Rd 



exp 



-R/Rd 



(6) 



2.2.2 Mass Distribution in Bulges 



The mass in bulges is spherically symmetric distributed ac- 
cording to an r 1//4 law: 



Sbuige = Sft,o exp 



where E(,,o gives the surface density of the bulge at the ef- 
fective radius, Rb, that contains half the projected light. 
However, the integral to turn this into mass as a function of 
radius has no analytic solution and so we use the following 
approximation instead: 



M bu i gc (< r) = Mi 



bulge 



rr / a 

Jo 



x(l + x) 3 



dx. 



(8) 



where Mb u i ge is the total stellar mass of the bulge and the 
auxiliary variable x = R/a where a = 0.56Rb- Integrating 
equation [8] the bulge mass inside a given radius r is given 
by 



M bu i gc (< r) = Mi 



bulge 



+ a 2 



(9) 



The pioneering work of lKormendvl l| 19771 ) demonstrated 
that there is a correlation between the effective radius (Rb) 
and the effective surf ace brightness (7) of elli pticals galax- 
ies. At the same time, iFaber fc Jackson! (|l976l ) showed that 
the luminosity of these galaxies is proportional to the ve- 
locity dispersion, with the derived relation L oc a 4 express- 
ing the intuitive notion that more luminous galaxies have 
higher velocity dispersions (and hence higher masses) . From 
these two correlations, and since L = 7riif (I), it follows that 
the effective radius Rb should be correlated with the veloc- 
ity dispersion a. Since the velocity dispersion is an intrinsic 
quantity of theoretical galaxies derived from the properties 
of their host halos, we can use it to determine the effective 
radius and the mass distribu tion of bulges. To do so, we use 
the proportionality law from iDiorgovski fc Davisl (|l987l ): 



10"' 
1CT 2 

o 
*E 

o 

Q. 

.? 10" 4 
& 
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M K - 5log 10 h 



-21 -20 



log(cr/kms _1 ) = 0.21 logiRb/h^pc) + 2.58. 



(10) 



Figure 1. Comparison of the predicted i<"-band luminosity func- 
tion at z = 0, from DLB07 (dashed red line) and the satellite 
disruption model (solid red line). The data poi nts represent the 
mean of a combined set of obser vations from (Col e et alj|2003 : 
iBell et al1 [2003; Jones ct al. 200^), with error bars reflecting the 
minimum and maximum estimates from the three data sets in 
each bin. 



3 PREDICTIONS FROM THE DLB07 MODEL 
WITH GALAXY DISRUPTION 

In order to better understand the changes introduced by 
the new physics, we start by plotting the new model pre- 
dictions leaving the parameters governing the other physics 
unchanged from the values in DLB07. In the next section we 
will use the techniques introduced in Paper 1 to find a new 
best fit, tuning the parameters governing the basic physics 
of the model in order to agree with a range of observations. 

3.1 The Luminosity Function 

In Fig. [T] we plot the K-band luminosity function for the 
new model with satellite disruption. The predictions at red- 
shift zero are compared with the values from the original 
DLB07 model and with a combination o f three different 
observational dat a sets l|Cole et al.l l200ll : IBell et all 120031 : 
I Jones et al.ll2006t ). respectively from 2DFGRS, 2MASS and 
6DFGRS. The final data points are given by the average of 
the maximum and minimum number density estimates in 
each magnitude bin, with errors o~i equal to half the differ- 
ence between them. The scatter in the combined observa- 
tional luminosity function represents the level of accuracy 
that we require from the model. 

The introduction of satellite disruption into the model 
has a considerable effect on the bright end side of the K- 
band luminosity function. In the framework of a hierarchi- 
cally growing Universe, these objects grow from material 
received during numerous mergers over their lifetime. With 
the inclusion of satellite disruption, a large amount of ma- 
terial which would otherwise end up in the brightest objects 
at redshift zero is transferred from satellites into the ICM. 
The result is an excessively low number density of bright 



© 2009 RAS, MNRAS 000.1111131 



Tidal Disruption of Satellite Galaxies 5 



1 .0 



0.6 



0.4 



0.2 



0.0 



• Best Fit with disruption 
o Best Fit without disruption 
■ Boldry et ol. (2004) 



9.5 



10.0 10.5 

log, (Mh 2 /M o ) 



1 1.5 



Figure 2. The predicted fraction of red galaxies as a function of 
stellar mass. The original lDe Luc ia & Blaizot ( 2007*) model (open 
red circles) is compared with the satellite di sruption mode l (filled 
red circles) and observational data from iBaldrv et al.l (|2004l ) 
(filled blue squares). 



objects in the K-bwad when compared to observations and 
the previous model. 

The number density of d warf galaxies is, i n its tu rn, less 
affected than we inferred in Irlenriques et al.l (|2008h . With 
the self-consistent implementation, disruption is no longer 
an instantaneous and dramatic process; instead galaxies lose 
material slowly as they follow their route into the central 
galaxy. Moreover, the merger time-scale is increased because 
satellites become smaller and less affected by dynamical fric- 
tion as they spiral inwards. Nevertheless, the decrease in the 
number density of dwarfs is still statistically significant, due 
to the large numbers of these objects contained in the Uni- 
verse. 

We emphasize that the poor fit results from the fact that 
we used parameters values tuned for the basic physics, with- 
out the inclusion of satellite disruption. This makes it easier 
to understand the changes introduced by the new recipe. In 
the next section, we will test if there is a region in param- 
eter space where disruption can produce an overall better 
agreement between the predicted and observed luminosity 
functions. 



3.2 Galaxy Colours 

In Fig. [2] we show the fraction of red galaxies as a function 
of stellar mass. We divide the gala xies into the two popula - 
tions using the selection criteria in Wcimnann et all l|2006h . 
(g - r) = 0.7 - 0.032 (M r - 5 log ft + 16.5), converted into 
a cut on the colour-stellar mass relation at redshift zero, 
(B — V) = 0.065 log(M*/i 2 /M ) + 0.09. The conversion 
from the g — r to t he B — V colour was done following 
iFukugita et all l|l996h . g-r = 1.05(5 - V) - 0.23. The frac- 
tion of red galaxies for differe nt mass bins is then compared 
with observations from IBaldrv et all l|2004h . The observa- 
tional masses based on the 'diet' Salpeter IMF (|Bell et al.l 



2003) were reduced by 0. 15 dex to agree with the IMF as- 
sumed in our SA model (jChabrierl l200St ) . The predictions 
from the model with satellite disruption are compared with 
the DLB07 model ( open red circles) and observational data 
from IBaldrv et aT] ((2004) (filled blue squares). The disrup- 
tion of stellar material from satellites significantly reduces 
the fraction of red galaxies over all mass ranges. 

This reduction comes about through the removal of pre- 
dominantly red stars from galaxies into the ICM. For central, 
mostly high-mass galaxies in halos, this reduces the mass of 
material accreted from merging dwarfs. 



3.3 The Black Hole-Bulge Mass Relation 

In Fig. [3] we show the black hole-bulge mass relation for 
the satellite disruption model. The colours follow the num- 
ber density of objects with blue representing low and green 
high density regions. In comparison with the predictions for 
the original model (Fig. 5 in Paper 1), there is an overall 
reduction in the masses of both bulges and of black holes. 

The reduction on black hole masses is determined by 
the black hole growth implementation in the model. The 
build up in black hole mass in the original DLB07 model is 
mostly due to the quasar mode: 



Am B HV 



/BH^sat/m-ccntral) "Icold 

l + (280kms- 1 /Kir) 2 



(11) 



This equation represents the fact that, during a merger 
event, the amount of cold gas driven into the central black 
hole depends on the instabilities created. The instabilities 
themselves depend on the mass ratio between the two merg- 
ing galaxies. Satellite disruption decreases the overall mass 
of satellites, reducing the instabilities created during merg- 
ers and hence the black hole growth due to cold gas accre- 
tion. 

Looking at Fig. [3] we see that the reduction in black 
hole growth is more significant than in bulge growth. The 
peak previously seen around a bulge mass of 10 10 Mq and a 
black hole mass of 10 s Mq has been reduced, with objects 
moving to lower black hole masses. 



3.4 Intra-Cluster Light 

On of the greatest advantages of introducing stellar disrup- 
tion into semi-analytic models is that it naturally explains 
the production of intra-cluster material. For the original set 
of parameters the disruption model predicts ICL fractions 
of approximately 22 per cent for the most massive clusters 
of galaxies M v ir = 10 15 Mq. In section [5771 we will compare 
these predictions with those of an optimized model for virial 
masses ranging from 10 12 to 10 15 Mq. 



3.5 The Metallicity of Stars 

Apart from the galaxy quantities that originally motivated 
the introduction of a new recipe into the SA model, there 
is another predicted property that is significantly improved 
by the disruption of satellite galaxies. In Fig. |4]we show the 
metallicity of stars as a function of the galaxy stellar mass. 
The new disruption model (solid red lines) is compared with 
the original DLB07 predictions (dashed red lines) and with 
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Figure 3. The black hole-bulge mass relation for the satellite 
disruption mod el (solid contours). T he red crosses represent ob- 
servations from lHaring fc Rixl ||2004| ) with the best fit to the data 
points given by the red line. The black lines represent the binning 
used to compare the model with observations in sections [4] and \5\ 



Figure 4. Comparison between the metall icity of stars in the 
satelli te disruption model (solid red lines), in lDe Lu ciafcBlaizot 
2007) (dashed red lines) and in observations from lGallazzT*et**aL 



20051 ) (blue squares and lines). For all the data sets, the central 



observations from iGallazzi et ail ((2005) (blue squares and 
lines) . The central line represents the median value of metal- 
licity in each mass bin (the blue squares for the observational 
data), while the upper and lower lines represent the 16th and 
86th percentiles of the distribution. From the figure it is clear 
that satellite disruption improves the agreement between SA 
predictions and observations, by increasing the metallicity 
of intermediate and high mass galaxies, M* > 10 10 M . The 
explanation for this effect is related to the decrease of the 
number density of the most massive galaxies. With the im- 
plementation of disruption, a considerable amount of ma- 
terial that would otherwise end up in these objects is now 
transferred into the ICL. In terms of the metallicity of stars, 
this results in massive objects receiving less low metallicity 
material from satellites, hence increasing their mean metal- 
licity. 



4 MCMC PARAMETER ESTIMATION 

In Paper 1 we have implemented MCMC sampling in the 
Munich SA model. This technique allows us to combine the 
constraining power from multiple observational data sets 
with a fast sampling of high-dimensional parameter spaces. 
By doing so we can verify the level of agreement with obser- 
vations and the relative weight of different observations in 
the final choice of the parameters in the best fit model, in a 
statistically-consistent way. Moreover, whenever reasonable 
agreement proves to be impossible, it helps us understanding 
whether there is a failure in determining the right parameter 
configuration, whether there is a fundamental problem with 
the underlying model, or whether the introduction of new 
physics is required. 

In the previous section we presented the predictions for 
the SA model with a self-consistent treatment of the dis- 



line represents the median value of metallicity in each mass bin 
(the blue squares for the observational data), while the upper 
and lower lines represent the 16th and 86th percentiles of the 
distribution. 



ruption of satellite galaxies. These were obtained using the 
parameters from the original model of DLB07, not adjusted 
to incorporate this new physical recipe. In this section, we 
will use the MCMC sampling techniques introduced in Pa- 
per 1 to search for the combination of parameters that gives 
the best fit for the satellite disruption model. We refer the 
reader to that paper for a full description of the statistical 
methods. The same MCMC implem entation is used, namely 
the Metropolis -Hastings algorithm l|Metropolis et af1ll953l ; 
lHastingall970r ) and a log-normal proposal distribution with 
a width that assures an overall acceptance rate in each chain 
between 10% and 40%. 

Due to the computational requirements, we will once 
more perform the sampling in a single volume of the mil- 
lennium simulation, representing 1/512 of the total simu- 
lation. Since we choose a representative file, with a lumi- 
nosity function and mean density similar to that of the to- 
tal volume, we are able to correctly constrain the proper- 
ties of galaxies in all mass ranges except for the most mas- 
sive objects (M* > 10 11 M ). The MCMC sampling is per- 
formed over ~30 000 steps and the output analyzed using 
GETDIST, which is pa rt of the COSMOMC software package 
IjLewis fc Bridldl2002l ). adapted to produce Id and 2d max- 
imum likelihood (profile) and MCMC marginalised (pos- 
terior) distributions. Once the best fit parameters are ob- 
tained, we then re-run the SA code over the entire Millen- 
nium to obtain the galaxies properties presented later in this 
paper. 

Using one file representative of the total volume sig- 
nificantly reduces the computational time required for our 
study. Nevertheless, the size of the calculations involved, 
even in that smaller volume, still make it a challenging task. 
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Table 1. SA model parameters from lDe Lucia &: Blaizod <|2007h ■ 
The first 6 parameters are frozen in our analysis at the values 
shown here. 



fb 


20 


Zi 


Tmcrgcr 


R 


Y 


0.17 


8 


7 


0.3 


0.43 


0.03 


asp 


*AGN 


/bh 


Edisk 


ehalo 


7ej 


0.03 


7.5 x 10~ 6 


0.03 


3.5 


0.35 


0.5 



In order to perform the sampling, the Cosmology Machine 
(COSMA) supercomputer supplied by Sun Microsystems 
was used, a machine based in the Durham University that 
is part of the Virgo Consortium facilities. 

4.1 Model Parameters 

As in Paper 1, in order to better understand the basic 
physics of the model, we choose to do our MCMC sampling 
only on 6 of 12 parameters in the model. On the top line 
of table [1] we show the frozen parameters corresponding to 
the baryon fraction (fb), the redshifts of beginning and end 
of reionization (respectively zo and z T ), the major to minor 
merger threshold (T merg e r ), the instantaneous recycled frac- 
tion (R) and the yield of metals (Y). On the bottom line, 
we show the values for the more fundamental parameters 
that we choose to sample. These are the star formation ef- 
ficiency («sf), the AGN radio mode efficiency (/cagn) , the 
black hole growth efficiency (/bh), the SN reheating and 
ejection efficiency (respectively edisk and thaio) and ejected 
gas reincorporation efficiency (7ej)- The physics governed by 
all the parameters are fully described in lCroton et al.l (|2006t ) 
as well as in Paper 1. 

4.2 Observational Constraints 

The comparison between model and observations will be 
done using the same data sets and statistical tests as before 
and since we have not introduced any additional parameter, 
the relative goodness of the original and the satellite dis- 
ruption model can be assessed directly. The ability or not 
for the MCMC to find a best fit model with a likelihood 
higher than before will tell us if the introduction of the new 
physical recipe is justified. 

The three observations used are the K-band luminosity 
function, the fraction of red galaxies versus stellar mass and 
the black hole-bulge mass relation. For a detailed descrip- 
tion of the observational data sets and the statistical tests 
employed, we refer the reader to Paper 1. 



5 RESULTS FROM THE OPTIMISED 
STRIPPING MODEL 

5.1 The Best Fit Semi- Analytic Model with 
Satellite Disruption 

We now present the results for the MCMC sampling when 
the SA model with satellite disruption is constrained by the 
three observational data sets combined. The likelihood for 
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Figure 5. Correlations between the 6 parameters analysed in the 
study for the SA model with satellite disruption constrained by 
three observational properties: the i<"-band luminosity function, 
the fraction of red galaxies, and the black hole-bulge mass rela- 
tion. The values of the parameters are plotted in log space, with 
the solid contours representing the 68% and 95% preferred regions 
from the MCMC (the posterior distribution) and the colours the 
maximum likelihood value sampled in each bin (the profile dis- 
tribution). The colour scale is normalized by the maximum like- 
lihood value of 0.15. 

the best fit model is given by the product of the likelihood 
from the three statistical tests: 

K — band) x £ ( c olour) X £(BH-Bulgc)- (12) 

Fig. [S] shows the allowed ranges and correlations be- 
tween the parameters sampled. As in Paper 1, the use of 
combined observations to produce one comprehensive data 
set in terms of galaxy formation properties restricts the pa- 
rameters to one small region with acceptable likelihood. 

The original model had a maximum likelihood value 
of just 0.037 when compared with the chosen set of obser- 
vational constraints, meaning it was formally ruled out at a 
2-cr level. However, the new satellite disruption recipe brings 
the model likelihood up to 0.15. Since we introduced no ad- 
ditional parameters, this increase in a factor of 4 in peak 
likelihood means that the inclusion of the new physical pro- 
cess seems to be favoured by data. 

In comparison to our analysis for the model without 
galaxy disruption, we see both an increase in the likelihood 
value of our best fit and in the allowed regions in parame- 
ters space. This is mainly caused by the changes that the 
new physical recipe produces on the black hole-bulge mass 



© 2009 RAS, MNRAS 000. 



8 Henriques et al. 



relation. The regions required in parameter space by the dif- 
ferent tests now have a much larger overlap. This means that 
the black hole and bulge mass build up in the new model 
is now consistent with the other galaxy properties analyzed, 
namely the if-band luminosity function and the fraction of 
red galaxies. 

In order to understand the changes introduced by satel- 
lite disruption we include two additional plots in our anal- 
ysis (the two top panels in Fig. [5} , the correlations between 
the AGN quasar mode parameter and the star formation 
efficiency, and the AGN radio mode and the disk reheating 
efficiency. 

An expected correlation, but one that was not seen in 
Paper f is shown in the upper left panel. A positive corre- 
lation is now evident between the quasar mode, responsible 
for the black hole growth, and the star formation efficiency, 
responsible for the bulge growth. Considering the physics 
governed by them, it comes at no surprise that an increase 
in one requires an increase in the other in order to maintain 
the black hole-bulge mass fraction of galaxies. The fact that 
this correlation was not present in the previous study is ex- 
plained by the considerably smaller regions with acceptable 
likelihood that we found before. 

For the same reason, the correlation between the super- 
nova reheating and ejection parameters, produced by the 
_K"-band luminosity constraint, is now clear in the lower left 
corner. As explained in Paper I this correlation keeps the 
virial velocity cutoff, above which SN feedback stops being 
effective, at a constant value: 

Vvir.O = ( — ) 2 ^SN- (13) 
\ Edisk / 

This means that this form of energy only suppresses star 
formation in the smallest objects where an excess compared 
to observations was previously seen. 

Accordingly, the lower and middle right panels show an 
identical behaviour for the parameters as in the T^-band 
sampling on Paper 1. An increase in the amount of gas 
ejected by SN needs to be balanced by an increase in the 
reincorporation time-scale and the star formation efficiency. 

5.2 Best Fit Parameters and Confidence Limits 

The best fit and confidence limits for the 6 free parameters in 
the model with satellite disruption, together with the pub- 
lished values from DLB07 and the best fit for the model 
without satellite disruption are shown in table [2] Despite a 
different model being used, the parameters from both the 
original DLB07 and the previous MCMC analysis remain 
close to the new 2a confidence limits. 

In comparison to the best fit obtained for the model 
without satellite disruption, the most dramatic effect in 
terms of preferred parameter values is the increase in the 
star formation efficiency. It represents a change in the frac- 
tion of cold gas transformed into stars in a disk dynamical 
time from 4% to 16%. Such an increase is required to repop- 
ulate the bright end side of the Tf-band luminosity function, 
largely affected by disruption. This effect is somewhat bal- 
anced by an increase in the mechanical heating from AGN, 
given by the product of &agn and /bh- 

The virial velocity cutoff, above which supernova stops 
being effective (Eq, ll3[) , is similar to that found for the model 



without disruption. This means that the new model still re- 
quires the supernova feedback to be effective only for small 
objects. However, the strength required from this form of 
feedback is now smaller (6.94 for the reheating and 0.33 for 
the ejection efficiency), which seems to be in better agree- 
ment with observations l|Martinll 19991 ). This means that dis- 
ruption does indeed help reducing the excess of dwarfs in the 
model and that an excessively strong feedback is no longer 
required. 

5.3 Predictions for the Best Fit Satellite 
Disruption Model 

In this section we analyze the predictions from the best 
fit model with satellite disruption. Contrary to the MCMC 
sampling, done in a single volume due to computational re- 
source limitations, the following predictions were obtained 
using the full Millennium Simulation. 

5.4 Galaxy Luminosity Functions 

In the left panel of Fig. [6] the A'-band luminosity functions 
from the best fit for DLB07 and the best fit for the satel- 
lite disruption models are plotted against the observational 
data set used to constrain the sampling. In comparison to 
the best fit obtained for a model without disruption the new 
model produces similar galaxy luminosities. However, for the 
model including satellite disruption, the good agreement at 
the low luminosity end is achieved with a less efficient super- 
nova feedback, which seems to be in better agreement with 
observations. 

In the right panel of Fig. [6] we show the predictions for 
the foj-band luminosity function. The best fit satellite dis- 
ruption model is compared with the best fit for DLB07 and 
two observational data sets. As for the A"-band, both best 
fit models with and without disruption show produce similar 
predictions, achieving an overall good agreement with obser- 
vations, except for the region around L*. We note however 
that the fej-band flux is highly dependent on the adopted 
dust model (which does not happen with the other proper- 
ties analyzed). For this reason, the excess just mentioned can 
be removed by adjusting the dust implementation, without 
affecting the agreement achieved for the other galaxy prop- 
erties. We choose not to do so, in order to make more clear 
the changes introduced by satellite disruption. 

5.5 Galaxy Colours 

Fig.Qshows the predictions for the galaxy colours in the best 
fit model with satellite disruption. The top panel gives the 
B-V colour-stellar mass relation (colour coded by the num- 
ber density of objects), while the bottom panel shows the 
fraction of red over the total number of galaxies as a func- 
tion of stellar mass. The introduction of disruption without 
changing the original DLB07 parameter values decreases the 
overall fraction of red galaxies, causing the model to under- 
predict the number of these objects except for L* galax- 
ies (Fig. [2}. As seen in the bottom panel of Fig. [7] , the 
MCMC optimization brings the model into agreement with 
observations for the entire plotted range. When compared 
to the best fit without disruption, the new model produces 
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Table 2. Statistics from the MCMC parameter estimation for the parame ters in the satellite disrup tion model. The best fit and 
marginalized confidence limits are compared with the published values from Dc Lucia fc Blaizoj |2007t ) and with the best fit values 
obtained without the inclusion of satellite disruption (Paper 1). 





DLB07 


Paper 1 


Disruption Model 


-2a 


-la 


+la 


+ 2(7 


<>SF 


0.03 


0.039 


0.17 


0.078 


0.13 


0.28 


0.53 


&AGN 


7.5 x lO" 6 


5.0 x 10~ 6 


5.3 X 10~ 6 


2.7 x 10~ 6 


3.7 x 10" 6 


6.2 x 10~ 6 


7.9 x 10~ 6 


/bh 


0.03 


0.032 


0.047 


0.030 


0.041 


0.061 


0.075 


Edisk 


3.5 


10.28 


6.86 


5.22 


6.33 


8.51 


10.11 


etalo 


0.35 


0.53 


0.33 


0.26 


0.31 


0.40 


0.46 


7cj 


0.5 


0.42 


0.13 


0.076 


0.12 


0.24 


0.30 




Figure 6. Comparison of the predicted A"-band (left panel) and bj-band (right panel) luminosity functions at z = from the best 
fit for DLB07 (dashed red line) and the best fit for the satellite disruption model (solid red line). On the left panel, the data points 
represent the observatio ns used to constrain the luminosities of galaxies in the MCMC parameter estimation llCole et al.ll2uQ]]:lBell et al.l 
l2003l; |j OIKS ct al. 2006). On the right panel, the ftj-b and luminosity functi on is compared w i th ob servations from 2DFGRS (green filled 
squares) and 6DFGS (blue open circles), respectively iNorberg et"al ] (|2002h and[j ones et al ] <|2006h . 



an overall better agreement with observations, reducing the 
systematic excess of small galaxies previously seen (open red 
circles) . 

More interesting is that the model now starts to produce 
an isolated population of massive red galaxies, as required by 
observations (top panel of Fig.[7]|. However, despite reducing 
the number of dwarf red galaxies, the disruption model still 
has an excessive fraction of these objects when compared to 
observations. 



5.6 The Black Hole-Bulge Mass Relation 

The changes introduced by satellite disruption in the pre- 
dicted black-hole bulge mass relation are the main reason 
for the better fit obtained in the new model. This means 
that the allowed likelihood region for the it-band luminos- 
ity function and for the fraction of red galaxies now produces 
black hole and bulge masses that agree with observational 
results. 

As explained in section 13.31 disruption has a large im- 



pact on both black hole and bulge masses. However, looking 
at the distribution of these masses in the best fit model, it 
becomes clear that the higher likelihood obtained is largely 
due to the reduction in black hole growth. The small sizes 
of the satellites will mean that smaller instabilities are cre- 
ated during mergers and that smaller amounts of cold gas 
are driven into the central black holes. 

In Fig. [8]we show the black hole-bulge mass relation for 
the satellite disruption model with the best fit parameters. 
The peak around bulge masses of 10 10 M and black hole 
masses of 10 s Mq has disappeared. This means that most of 
the intermediate mass bulges and black holes are now on or 
below the observational red line, where most observational 
points are located. The binomial test used produces a prob- 
ability three times higher than before. 

5.7 The Intra-Cluster Light 

In Fig. [9] we plot the predicted ICL from our best fit model 
with satellite disruption. The solid line represents the me- 



© 2009 RAS, MNRAS 000,[THT3l 



10 Henriques et al. 




10 



10.5 11.0 
log, (Mh 2 /M s ) 



12.0 



Horing & Rix (2004) 



12 



13 



lo g( m Bukje/ M s) 



1.01 — r 



0.4 



0.2 



0.0 



9.5 



T 



i 



• Best Fit with disruption 
o Best Fit without disruption 
■ Boldry et al. (2004) 



10.0 10.5 

log, (Mh 2 /M o ) 



1 1.0 



1 1.5 



Figure 7. The top panel shows the B-V colour-stellar mass re- 
lation for the galaxies in the best fit satellite disruption model. 
The solid line represents the division between the red and blue 
populations in Wein mann et ahl (2006). The predicted fraction of 
red galaxies as a function of stellar mass is showed in the bot- 
tom panel. The best fit for the original DLB07 model (open red 
circles) is compared with the best fit for the satelli te disruption 
model (filled red circles) and observational data from Bal drv et al.l 
\2004 ) (filled blue squares). 



dian of the MiCL/Mtotal distribution in each bin while the 
dashed line gives the same relation for the disruption model 
with the original set of parameters. 

The predicted ICL fraction for groups with M v i r > 
10 13 M Q has a mean of ~ 18%. This value is compatible with 
observations that detect ICL fractions between 10% and 40% 
dBernstein et al . 1995; G onzalez et al.|[2 000: Fcld meier et al.l 
120021 , 12004 iGonzalez et al.ll2005l ; IZibetti et al.ll2005T l. The 
prediction decreases systematically has we move to smaller 
virial masses reaching a mean of « 7% for M v i r = 10 12 Mq . 
The difficulty in distinguishing between ICL and that from 



Figure 8. The black hole-bulge mass relation for the best fit 
satellite disruption mo del (solid contours). T he red crosses repre- 
sent observations from lHarine fc R ix ( 2004) with the best fit to 
the data points given by the red line. The black lines represent 
the binning used to compare the model with observations. 
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Figure 9. Fraction of the mass in the ICM over the total stellar 
mass of the group as a function of virial mass. The blue dots are a 
representative sample of the total galaxy population in the best fit 
model with satellite disruption. The solid and dashed lines repre- 
sent the median of the MiciV-Mtotal distribution for the satellite 
disruption model with the best fit and original parameters respec- 
tively. 



the halo of the central galaxy (the former can be regarded as 
an extension of the latter) , means that observational data is 
only available for large groups and has a considerable scat- 
ter. 
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Figure 10. Comparison between the metallicity stars in the 
best fit for the satellite disruption model (solid red lines), in the 
best fit for DLB07 ( dashed red lines) and in observations from 
Gallaz zi et al ] lj2005l l (blue squares and lines). For all the data 
sets, the central line represents the median value of metallicity in 
each mass bin (the blue squares for the observational data) , while 
the upper and lower lines represent the 16th and 86th percentiles 
of the distribution. 

5.8 The Metallicity of Stars 

In Fig. [10] we show the metallicity of stars as a function of 
the galaxy stellar mass. We remind the reader that this was 
not one of the observations used to constrain the model. The 
best fit disruption model (solid red lines) is compared with 
the best fit for the ori ginal DLB07 (dashed red lines) and 
with observations from iGallazzi et all (|2005T ) (blue squares 
and lines) . As in Fig. [4j the central line represents the me- 
dian value of metallicity in each mass bin (the blue squares 
for the observational data), while the upper and lower lines 
represent the 16th and 86th percentiles of the distribution. 

As for the best fit model without disruption, the in- 
crease in supernova feedback required by the MCMC lowers 
the metallicity of stars for the low mass galaxies in the best 
fit for the disruption model. In comparison to the original 
model (Fig. [3), for this mass range, the higher SN feedback 
from both best fit models increases the fraction of metals 
in the gas phase, consequentially lowering the fraction of 
metals in stars. 

As explained in section 13.51 the introduction of disrup- 
tion increases the metallicity of massive objects, by reducing 
the amount of low metallicity material that they receive from 
mergers. This means that the distribution of metals in the 
model is now in close agreement with observations over the 
entire mass range plotted. 



6 DISCUSSION 

In Paper 1, we introduced a new approach to galaxy forma- 
tion modelling. Combining the MCMC sampling techniques 
with SA models allowed us to gain insight on the physical 



importance of the different galaxy parameters. With the in- 
troduction of a new physical recipe, presented here, both 
the normalisation and shapes of the acceptable likelihood 
regions contain useful information that can be used to dis- 
criminate between models. 

In this work we implement a new physical ingredient 
in the SA recipe: the stripping of satellite galaxies by tidal 
forces during merging events. The co ncept had already been 
introduced in lHenriques et all l|2008l ). but only as an a pos- 
teriori study. That approach, despite enabling us to gain 
insight on the impact that this process would have in the 
population of dwarfs, did not allow us to study the effect 
of stripping in a self-consistent way, and in particular the 
effect on the most massive galaxies. 

The self-consistent implementation that we describe in 
this chapter makes it possible to study the impact of dis- 
ruption on the properties of galaxies of all types. Namely, it 
allows us to study the loss in mass by satellite galaxies, the 
slower build up of central galaxies (since they receive less 
material from satellites) and the growth of the ICM com- 
ponent, previously neglected. On top of that, the MCMC 
sampling allows us to learn if the new process is favoured or 
not by observations. 

In comparison to the original model, the best fit likeli- 
hood of the model with disruption is four times higher with 
respect to the K-band luminosity Function, the fraction of 
red galaxies and the black hole-bulge mass relation. 

• Since we did not introduce any additional parameter 
to model disruption, this means that the inclusion of tidal 
disruption of stellar material from satellite galaxies during 
mergers seems to be favoured by observations. 

• Moreover, it means that the new model is now formally 
consistent with the combined observational data set we used. 

• The higher likelihood value and the larger allowed like- 
lihood regions in parameter space of the new best fit model 
are mainly determined by the changes produced by satellite 
disruption on the black hole-bulge mass relation. This means 
that the growth of bulges and black holes is now consistent 
with the properties of the galaxies as a whole. 

• The new best fit has a considerably higher star forma- 
tion efficiency in order to correct the reduction in the num- 
ber density of massive objects caused by disruption. This 
is balanced in massive galaxies by an increase in the AGN 
feedback efficiency. At low masses, a less effective supernova 
feedb ack is now req uired, in better agreement with observa- 
tions i|Martidll999f ). 

• Although they were not used as constraints, the MCMC 
sampling for the new model kept the agreement with the 
observational foj-band luminosity function and significantly 
improved the overall shape of the metallicity distribution of 
stars. 

• Finally, the introduction of disruption allows us to fol- 
low the build up of the ICM, which represents about 18% 
of the total light in clusters, in agreement with observations 
(between 10% and 40%). 

6.1 Future Challenges 

Despite the higher likelihood found for the best fit model, 
SA predictions are still far from exactly reproducing obser- 
vations. From our analysis, the biggest challenge for the SA 
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model is still to reproduce the properties of dwarf galaxies. 
With the introduction of satellite disruption, the model is 
able to reproduce the number density of these objects with 
a less efficient supernova feedbac k, which is in better agree- 
ment with observational studies (|Martinlll999l) . 

Nevertheless, this type of galaxies remain predom- 
inantly red, in disagreement with observational studies. 
Stripping of stars does help reducing the number of red 
dwarfs, but is only part of the solution to the problem. An- 
other important factor might be the treat ment of the i mpact 
of mergers on the gas phase. For example. iFont et alj (I2008T ) 
consider a model in which the gas of satellite galaxies is con- 
tinually, instead of instantaneously stripped, enabling them 
to form stars and remain on the blue sequence for longer. 

The higher likelihood found for the best fit model with 
satellite disruption seems to be directly related to a bet- 
ter agreement between the black hole-bulge mass and the 
other two constraints. This fact is apparently related to the 
slower build up of black holes relative to bulges and the 
overall slower build up of both components. However, we 
emphasize the points raised in Paper 1 about the simplic- 
ity of the black hole growth model, which still appears to 
neglect important features such as the impact of the quasar 
mode feedback on galaxy properties. In addition, from an 
observational point of view, larger and more robust data 
sets of black hole properties are needed. 

The predictions obtained for the ICL are in good agree- 
ment with observations. However, the available data is still 
subjected to large error bars, due to the difficulties in distin- 
guish between the light from the brightest cluster galaxies 
and the ICM itself. Also, this component is still only de- 
tectable in large groups and clusters due to its low surface 
brightness. If future improvements in observations make it 
possible, this property could be used to directly constrain 
the model with satellite disruption. 
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